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Abstract 

This study aimed to elaborate and develop the existing methods used in inventorying lowland forest and adapt them 
for mountainous forests. 

The inventory of mountain forests with the use of CIR aerial photographs in this study relied on the imple- 
mentation of the 3D (three-dimensional) methods and it was characterised by two phases. In the 1 st phase the 
auxiliary variables (e.g. crown cover, height of dominant trees, density) were measured on the aerial photographs 
(355 circular plots). Variable of interest (standing volume) was recorded during the terrestrial survey on the cor- 
responding ground sample plots. A statistical relationship was established between the variable of interest and the 
auxiliary variables by means of multiple regressions. In the 2 nd phase selected auxiliary variables were measured 
on the enlarged set of 2772 plots on the aerial photographs only. The variable of interest (standing volume) was 
computed for the whole forest area by using the regression model developed in the 1 st phase and the aid of the 
variables measured in the 2 nd phase. 
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Introduction 

Forest inventory as the key planning and controlling 
tool in forestry is based on the data collection, which 
can be separated in 1) (measured) primary data, 2) (cal- 
culated) secondary data and 3) estimated data and re- 
sults in information about the forest stands attributes 
such as: tree species, tree species composition, breast 
height diameter (DBH), tree height, number of trees, 
basal area, age and more importantly the estimation of 
standing volume and volume increment. The process 
of forest inventory which is carried out to provide such 
data mainly is conducted by using traditional methods, 



which means all the necessary primary data are meas- 
ured on the ground (Akca et al. 1996). 

Traditional way of forest inventories appears to be 
quite time consuming, plus requires a huge amount of 
human labour and culminates in high costs. Therefore 
the implementation of new technologies combined with 
the traditional way of forest inventory could be more 
promising in reducing the time for accomplishing such 
tasks and simultaneously the costs, thus increasing the 
efficiency of forest management. 

The inventory of mountain forests with the use of 
aerial photographs relies on the implementation the 3D 
(three-dimensional) methods. This is much more rele- 
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vant in mountain forests where often aerial photographs 
have been used to reduce the effort of data collection. By 
increasing the technical standards, especially in com- 
puter science during the last decades, 3D methods can 
be implemented in such processes, however this does 
not substitute completely the traditional ground-based 
forest inventories. There are needed ground-based in- 
ventory data sets to compare with the aerial data in 
order to determine the accuracy e.g. standing volume 
estimation (Kohl et al. 2006; Michael et al. 2008). 

The aerial photographs were recorded on colour 
infrared (CIR) composition. Such a composition is cho- 
sen based on the fact that the photosynthetic activity of 
the vegetation has its maximum radiation reflectance in 
the near infrared from 700 to 950 nm (Lillesand et al. 
2004). 

CIR plays a key role in discriminating the tree spe- 
cies on the aerial photographs because it provides a di- 
versity of multi-spectral information (Ross 1970). Be- 
cause of the variations in the photosynthesis patterns, 
the radiation reflectance on the aerial photograph tends 
to be different on various types of growing vegetation 
thus enabling us to distinguish between tree species 
(Ross 1970). 

Remote sensing technology and in particular 
the use of aerial photography, in the last decades has 
proven to be very effective in forest management ap- 
plications such as forest inventory or forest monitoring. 
Furthermore it has significantly enhanced our ability to 
improve the inventorying of mountain forests. A great 
deal of credit can be ascribed to the incorporation of 
the DTM (Digital Terrain Model) with the aerial photo- 
graphs which has a great influence on the accuracy of 
tree height measurements and subsequently on the over- 
all volume estimation results (Baltsavias 1999; Hyyppa 
and Hyyppa 1999). 

The aim of this study was: 

- to elaborate and develop existing methods used in 
inventorying lowland forests and adapt them for 
mountainous forests, 

- to increase efficiency and precision to provide reli- 
able forest inventory data which are comparable and 
consistent, 

- to further develop methods and analysis techniques 
to facilitate the use of aerial photographs for esti- 
mating the biophysical forest attributes. 



Material and methods 
M ountainous forests 

Around 20% of all forest in Europe are mountain for- 
ests and there is an increase awareness of the manage- 
ment and assessment of these forest ecosystems (EOMF 
2000 ). 

Naturally, mountainous forests are associated with 
topographic variations and fluctuations resulting in 
unique landscapes. Generally the occurrence of forest 
communities changes depending on the elevation, slope 
and aspect consequently stimulating the prevalence of 
certain species (Gluck 2001). 

Elevation which refers to the altitude of the stand 
above sea level, determines in particular the vertical 
distribution of the tree species. Slope is another attrib- 
ute of mountainous forest which is strictly related to 
the variations of the terrain (e.g. hill top, hillside, plain, 
valley, ridge, etc). Aspect on the other hand refers to 
the geographical orientation of the hill-slopes, in other 
words the direction which the forest stand is facing (e.g. 
east, west) and affects the horizontal tree species distri- 
bution (Jin et al. 2008). 

The study area was the Stolowe Mountains Nation- 
al Park (Stolowe Mts NP) which is located in south-west 
of Poland. Stolowel Mountains National Park is com- 
posed of different forest stands, where spruce stands 
dominate among the others. From the total forest area of 
5,671.8 hectares Norway spruce ( Picea abies L.) stands 
represent more than 80% (J^dryszczak and Miscicki 
2001). There are other tree species considered as the 
ruling species like Scots pine ( Pinus sylvestris L.), Eu- 
ropean larch ( Larix decidua Mill.), silver fir ( Abies alba 
Mill.), European beech ( Fagus sylvatica L.), sycamore 
(Acer pseudoplatanus L.) and birch (Betula pendula 
Roth, and B. pubescens Ehrh.). However often the area 
which they occupy is rather low, because the only tree 
species whose share exceeds 5% except spruce is beech, 
the other tree species share vary from 0.1 to 5%. 

Two-phased method, as the term implies is char- 
acterised by two phases. In the 1 st phase data are col- 
lected from two samples of the same population. In the 
1 st sample which is represented by the aerial plots, the 



1 “Stolowe” in Polish means “table-like” therefore Stolowe 
Mountains stands for Table Mountains referring to the shape of these 
mountains 
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auxiliary variables (e.g. crown cover, tree height, den- 
sity) are measured whereas in the 2 nd sample the varia- 
bles of interest (e.g. DBH, tree height, age) are recorded 
during the terrestrial survey. The intent is to establish 
a statistical relationship between the auxiliary variables 
and the variables of interest by means of linear regres- 
sion (Kohl et al. 2006). In the 2 nd phase, using auxiliary 
variables which were chosen based on the calculations 
carried out in the 1 st phase, an extended population of 
sample plots are measured. Afterwards variables of in- 
terest (e.g. standing volume) were computed by exploit- 
ing the regression estimators derived from the 1 st phase 
and results of the 2 nd phase. 

The mean value is calculated for the inventory unit 
according to formula given by (Loetsch and Haller 1964). 

v=v N +b m {x 1 -x 2 ) (1) 



R 2 - coefficient of determination, 
p - number of used variables, 

N - population size, 

- size of the aerial sample plot (Phase I), 

- size of the ground sample plot (Phase II). 

Measurement of ground sample plots 

The terrestrial measurements were performed during 
July 2009 by the team of TAXUS SI. Measurements 
on the terrain consisted on a set of 355 circular perma- 
nent sample plots with a radius of 12.62 m. Permanent 
circular sample plots were distributed in a systematic 
grid arrangement each 400 meters on the entire area of 
the Stolowe Mts NP. From these plots the ground refer- 
ence data were collected and afterwards the necessary 
analysis and calculations were undertaken to estimate 
the ground volume. 



where: 

v., - mean value calculated from terrestrial measure- 
ments in phase I, 

b y( - the slope coefficient between v (measured on the 
ground plots) and x (measured on the aerial plots) 
computed from the phase I, 

Tj - mean value calculated from all aerial sample plots 
(phase I + II), 

x, - mean value calculated from the aerial sample plots 
of phase I. 

Slope coefficient is calculated: 




where: 

r - correlation coefficient, 

S - covariance characteristics of v andy, 

S - standard deviation of x. 

X 

In order to determine the square of the standard er- 
ror of mean value S 2 following equation of Khan and 
Tripathi (1967) was used: 
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where: 

S 2 - variance of the variable of interest in the first phase, 



Acquisition of aerial photographs 

Before the acquisition of the images one of the most im- 
portant steps is to determine the image scale. The scale 
depends on the objective and the accuracy require- 
ments. In this study photogrammetric data were com- 
posed of 12 strips of aerial photo stereo pairs acquired 
on 03/07/2008 by SWISS Photo Company. 

Digital images referred as photographs were ac- 
quired at a scale of 1:10000. The images were captured 
by airborne with a Vexcel UltraCam D camera using 
objective with a focal length 101.4 mm which registered 
four channels: B, G, R and 1R. In this study CIR com- 
position was used and spatial resolution of the images 
amounted a nominal pixel size of 13 cm. The aerial pho- 
tographs were taken at an altitude of 2100 metres above 
sea level with a side lap 60% and end lap 30%. 

Preparation of aerial imagery 
for photogrammetric measurements 

The interpretation of the aerial photography is based on 
the photogrammetric procedures which are executed by 
using “DEPHOS” digital photogrammetric software. 
During the process of the analysis and observation pairs 
of aerial photographs with stereoscopic overlap areas, 
which can be perceived in a three dimensional format, 
have been analysed and interpreted based on a scale of 
1:10000. In this study the used coordinate system was 
PUWG 1992. 
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Measurements of aerial sample plots 

Tree species identification and composition 

Different tree species are associated with various inten- 
sities of reflection which has an influence on the colour, 
tone and texture of the aerial photographs. This differ- 
ence allows to distinguish between tree species, how- 
ever the difference among coniferous species or among 
deciduous is more subtle than in case of coniferous ver- 
sus deciduous. 

It is important for interpretation to achieve CIR 
photographs as close as possible to the typical one and 
to amplify features of observed objects. For this reason 
LUT (Look up Table) in DEPHOS was used. 

After identifying all tree species which compose 
the plot, the second task is to estimate tree species com- 
position. Tree species composition is expressed in per- 
centage and represents the share within the plot which 
ranges from 1-100%. It should be noted, however, that 
the estimation could vary among different observers. 
Alternating between the observation of aerial photo- 
graphs and forest stand maps of the study area has im- 
proved tree species classification. 

Tree height measurement 

Tree height is one of the key parameters in forest invento- 
ry and can accurately be estimated from the aerial photo- 
graphs (Duvenhorst 1995). Tree height measurements on 
the ground are considered to be more difficult and time 
consuming than measuring the DBH, especially in dense 
deciduous stands where sometimes it is impossible to see 
the top of the tree. Therefore retrieving tree heights from 
aerial measurements is of great advantage in accelerating 
the process of forest inventory. 

Tree height is represented or is obtained by the dif- 
ference of tree top and the terrain elevation at the base 
of the tree. The procedure goes as follows: once the tree 
of interest is located by placing the 3D cursor on the 
assumed tree apex (z-tree-top) its position is recorded 
(x,y,z coordinates). The coordinates at the base of the 
tree correspond with the ground elevation (zground) 
from the DTM (Digital Terrain Model). Therefore sub- 
sequently the height of the tree is attained as the result 
of the subtraction of ground elevation from the (z) coor- 
dinate of the tree apex. 

The average plot height was obtained from the 
share of average height of each tree species within the 
plot (eq. 4). 



Average plot height 

Share of each 
tree species 

in % x average height ^ 

_ of each tree species 
vg - eig 100 

Height reduction was applied in order to avoid 
the estimation of the non merchantable volume. Trees 
as high as 7 meters do not produce assortments which 
could be considered merchantable therefore such trees 
are not inventoried during the field survey. 

Density measurements 

Average density is one of the parameters which is de- 
rived from the measurements on the aerial photographs 
and is one of the variables which might have been used 
for the further calculations of the standing volume. On 
each plot are measured 6 distances between the six clos- 
est neighbours from the tree which is located closest to 
the centre plot. In case when there are not enough trees 
within the plot the measurements can take place outside 
the plot. 

The formula used to calculate the average density is 
displayed below: 

Average Density 

12 200 

d= 7 V ( 5 ) 

A + ^ + - + A, 

V 6 

where: 

D - average density, 

L.+L.+ ...+L. - distances between the trees. 

12 o 

Density is expressed in trees ha’ 1 , and by multiply- 
ing the average height with the density another param- 
eter was obtain such us “Sum of heights”. 

Sum of heights 

Sum of heights = avgheight x density (6) 

where: 

avgheight - average height at plot level, 
density - average density at plot level. 

Sum of heights is a variable used in building up the 
model for the estimation of the standing volume. 
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Crown cover estimation 

The crown cover of a sample plot is determined as the 
percentage of the area which is covered by the crowns 
vertical projection of the dominant and co-dominant 
trees in a given plot. 

Basically this can be estimated through either 2D or 
3D observation and the estimation might differ slightly 
from one observer to another. The percentage of the 
crown coverage is dependent on the tree density, age, 
tree species (in case of shade tolerant species it might 
be more dense) and tree species composition. Crown 
cover is an important feature in describing the horizon- 
tal stand structure. 

Green volume calculation 

“Green volume” was calculated for each tree species 
and was derived as a sum of multiplication between 
average plot height, crown cover and the share of the 
particular tree species. 

Calculation of green volume for particular species 
GV = avg _height x CC% x tree _species% (7) 

where: 

avg_height - plot average height, 

CC% - plot crown coverage, 

tree_species% - share of particular tree species within 
the plot. 

Total green volume for the plot was calculated as 
a sum of the volume of each tree species present in the 
sample plot. Green volume at this stage represents the 
imaginary cylinder with a base diameter equal to the 
tree trunk diameter, therefore it has to be corrected. 

Age estimation 

Age is a very informative attribute about the forest, and 
it can be either taken from the available forest manage- 
ment reports or determined through counting the an- 
nual rings of the trees. In case of aerial sample plots it 
was checked whether age given in the stands descrip- 
tion corresponded with the height of trees and its crown 
diameter. 



Ground volume calculation 

Standing volume is one of the most important attributes 
of the forest and it is a parameter of a great interest from 
the forest management and planning point of view. 

Tree volume was a function tree DBH, tree height 
and volume coefficients of a particular tree species. 
Ground volume was calculated based on the measure- 
ments from the field investigation and was calculated 
for each inventoried individual tree. Plot volume was 
derived as a sum of the individual trees within the plot 
and extrapolated for the extended area. The volume for 
each tree species was calculated by using the formula 
Schumacher-Hall (Loetsch et al. 1973): 

Ground volume calculation 

Y — JQO+Mogrf u +Clog* 

where: 

a,b,c- coefficients which are depended on the tree 
species and the age group, 
c/j - DBH (diameter at breast height), 

h - tree height. 

Coefficients a, b and c used in Schumacher-Halfs 
volume equation for each tree species and correspond- 
ing age groups were based on the yield tables (Grudner 
and Schwappach 1952). 

Results obtained from the ground measurements 
were used as reference data and also as construction 
part in the development of regression models, particu- 
larly ground volume. 

Calculation of variables on aerial sample plots 

During the aerial measurements three sets of variables 
were computed. The first set of variables was derived 
from the direct measurements e.g. height, crown cov- 
er, average distance between trees, elevation. Belong- 
ing to the same set additional variables were comput- 
ed e.g. average height, sum of heights, trees density, 
green volume. Second set of variables was derived 
from the multiplication of the variables of the 1 st set, 
all variables were squared up producing new variables 
(e.g. CC% 2 , avg height 2 , Greenvol 2 ). Third set of the 
variables was a combination of first and second vari- 
able set e.g. (CC% x Av_Height), (Elevation x Height), 
(Green vol 2 x Density). 



Folia Forestalia Polonico, series A, 20 12, Vo I. 54 (2), 123-133 



128 



Besim Hoxha 



Calculation of the relationship between ground 
and aerial variables 

By using correlation analysis we observe what the asso- 
ciation between several variables individually and col- 
lectively is and later use them in multiple regressions. 
All correlation are described in a 14 x 14 correlation 
matrix (Tab. 1). 

What can be pointed out is a high positive relation 
between AH and Green Vol (0.85) and AH and Age 
(0.60). Furthermore it could be mention that a negative 
relationship exists between AH and Elevation (-0.56), 
due to the fact that species show less height growth 
by rising up the altitude. Since all the variables are 
considered statistically significant they can be used in 
building up the regression model in order to derive es- 
timators which will be used in predicting the standing 
volume. 

Multiple regression models were built to gener- 
ate equations in order to relate the parameters derived 
from the measurements on aerial sample plots (e.g. tree 
height, crown cover, density, green volume and eleva- 
tion) and volume derived from the ground plots. 

Multiple regression analysis was used to find the 
best fitting model for predicting standing volume. Dif- 
ferent models were tested in order to determine param- 
eter’s values of various functions based on the data set 
where ground volume was the dependent variable. 



Measurement of the extended part of aerial 
sample plots 

In the 2 nd phase measurements were carried out on 2772 
aerial plots distributed on the entire area of Stolowe 
Mountains National Park. During this phase only a few 
selected parameters were measured, unlike the 1 st phase. 
The selection of the parameters to be measured was de- 
termined based on the regression model which was used 
to estimate the volume in the 1 st phase. The variables of 
interest which were measured on the aerial photographs 
in the 2 nd phase were selected based on their influence 
in final regression model. The auxiliary variables with 
no strong influence on the statistical relationship were 
excluded and not measured in the 2 nd phase. 

Volume computations for the plots measured in 
the 2 nd phase were done by using the regression mod- 
el which was developed from the combination of the 
auxiliary variables and variables of interest in the 1st 
phase. 

Variants of inventory with the use of aerial 
photographs 

Once the volume was computed for each plot of the en- 
tire population in the Stolowe Mts NP, these plots were 
stratified into age classes. Computations of standing 
volume were carried out for each age class by using 
three different inventory approaches: 

- aerial method (two-phased method), 

- classical ground method and 

- combined method. 



Tab. 1 . Correlation coefficients between the main variables which is described by the 14114 correlation matrix. With: AH - Aver- 
age height; Elevation - Elevation of the sample plots; Age - Species-specific average age; Green_Vol - species-specific predicted 
green volume; N - number of plots. **. Correlation is significant at the 0.01 level (2-tailed) 



Variables 


Parameters 


Variables 


AH 


Elevation 


Age 


Green Vol 


AH 


Pearson Correlation 


1 


-0.506** 


0.600** 


0.850** 


Sig. (2-tailed) 




0.000 


0.000 


0.000 


N 


355 


355 


355 


355 


Elevation 


Pearson Correlation 


-0.506** 


1 


-0.250** 


-0.527** 


Sig. (2-tailed) 


0.000 




0.000 


0.000 


Age 


Pearson Correlation 


0.600** 


-0.250** 


1 


0.461** 


Sig. (2-tailed) 


0.000 


0.000 




0.000 


Green Vol 


Pearson Correlation 


0.850** 


-0.527** 


0.461** 


1 


Sig. (2-tailed) 


0.000 


0.000 


0.000 
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When using aerial inventory method, the volume 
for each age class is computed based on the same meth- 
odology as for the entire population (two-phased meth- 
od), taking into account the surface unit. Accordingly 
computations and statistical tests were performed to de- 
termine mean volume, coefficients (R, R 2 , b), standard 
and relative errors. 

Before applying of combined method, an analysis of 
covariance was used to compare the slope coefficient (b) 
between particular age classes. The corresponding co- 
efficients for each age class are switched for a common 
coefficient derived from the whole population (all age 
classes together). This is done in order to avoid the big 
differences among age classes within the allowed range. 
Practically the critical value of F-coefficient should not 
exceed 1.55 and in our case it was 1.066 therefore we 
were allowed to apply a common coefficient b = 1.029 
which represented the slope degree. Note, when using 
combined method in volume computation, volume is es- 
timated with some error, therefore can be used only for 
comparison not as a standard method. 

In the classical inventory method, for each individual 
tree species volume is calculated based on the formula 
of Schumacher-Flall (see Equation 9) and afterwards 
summed up at a plot level. Stratification into age classes 
has been applied only for permanent ground plots. 

Statistical tests at P = 0.05 significance level were 
run to verify and examine the accuracy of the three 
methods in order to determine their performance. 

Results 

Relationship between ground 
and aerial variables 

The selected model is a combination of the independent 
variables such as logarithmic Average height, logarith- 
mic Green volume, and the square of the Green volume 
with Ground volume as dependent variable. Logarith- 
mic method was used in order to avoid negative values 
because it would not make sense to have negative tree 
standing volume. 

Volume calculation (Model 3) 

Y = — 1.412 + 0.220 xH//_ log + (9) 

+ 2.108 x Green _ log -0.258 x Green _V 2 



where: 

Y - estimated volume, 

AH _ log - logarithmic Average height, 

Green _ log - logarithmic Green volume, 

Green_V~ - the square of the Green volum. 

Regression models were compared to each-other 
for their statistical parameters (e.g. R 2 , R 2 dj , residual 
of mean squares) and their simplicity. There are a few 
criteria which are used to compare regression models 
found in literature. They are: coefficient of determina- 
tion R 2 , the larger the value the stronger is the relation- 
ship. Adjusted determination coefficient (R 2 dj ) is more 
comparable if the models are built by different numbers 
of variables. Residual of mean squares (MSE) is another 
criterion, a model with smaller MSE is more preferable. 
However R 2 d . and MSE quite often infer the same contu- 
sion (Rawling 1988). 

Tab. 2 summarizes the results derived from 3 re- 
gression models. Based on the criteria mentioned above, 
3 rd model was selected to predict the volume because it 
provided almost the same accuracy and it was built with 
less variables than the rest of the models. The idea be- 
hind is to reduce the time which is required to measure 
these variables on the aerial photographs. The differ- 
ence in accuracy between the 3 rd model and the other 
two was negligable. 

Tab. 2. Comparison of the best regression models selected 
for the prediction of standing volume. Models were examined 
based on their statistical parameters e.g. determination 
coefficient (R 2 ), adjusted detennination coefficient (R 2 d ), 
residual of mean squares (MSE) at a signficance level * 



Model 


Auxiliary 

Variables 


R 2 


R 2 

adj 


MSE 


Signifi- 

cance 


1 


Green Vol; 
AgeGreenV; 
AH Green; 
Elev_AH 


0.747 


0.745 


11.232 


<0.001 


2 


Green Vol; 
Age_GreenV; 
AH Green; 


0.746 


0.744 


11.273 


<0.001 


3 


AH Jog; 

Greenjog; 

Green_V 2 


0.726 


0.724 


11.184 


<0.001 



The use of the regression model built from the com- 
bination of average height, green volume and squared 
green volume revealed a high correlation and a linear 
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realationship. The coefficient of determination amounts relation R = 0.852 indicating a quite high relationship 

to R 2 = 0.729, whereas the coefficient of multiple cor- and efficient variable estimator selection. Results of the 




£ 
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Age Classes 

Fig. 2. Average standing volume for all age classes in forested area of Stolowe Mts NP, calculated by two-phased, ground and 
combined method (bars represent confidence limits Cl at P = 0.05) 
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correlation between (ground volume) and (predicted 
volume) derived from the aforementioned model are 
plotted in Fig. 1. 

Standing volume in age classes 
and National Park 

Standing volume was computed for each age class by 
using three methods and a comparison of the results is 
illustrated in the Fig. 2. 

By the examination of the results (Fig. 2) it was ob- 
served that there was an underestimation of the volume 
from the aerial photographs for the old stands over 90 
years, especially stands over 141 which show the high- 
est difference. 

The average volume for all age classes calculated 
by contemporary ground method produced the follow- 
ing results 421 ± 15 m 3 ha' 1 and a relative error 3.6% at 
P = 0.05 statistical significant level whereas two-phased 
method was quantified to 403 ± 11 m 3 ha' 1 (i.e. 2.8%). 
Combined method on the other hand 370 ± 13 m 3 ha 1 
(i.e. 3.5%). We can see that generally there is an un- 
derestimation of the volume when using the combined 
method especially in comparison with ground method. 



The average volume for all stands without dividing 
them into age classes calculated by two-phased method 
was 339 ± 13 m 3 ha' 1 (i.e. 3.9%). 

Tab. 3 gives insights about the age classes with the 
corresponding mean volume, number of plots measured 
on the aerial photographs, on the ground and for the en- 
tire data set (3127 plots). 

The correlation coefficient R for stands from 1-60 
years old ranges not in a specific order from 0.736 to 
0.967, whereas for middle aged stands 60-100 the 
value of R coefficient ranges from 0.628 to 0.913, and 
for stands over 100 years old the achieved R value was 
rather higher than the previous ones ranging from 0.773 
to 0.913. Flowever, the mean correlation coefficient R 
had a value 0.852 which indicates a high relationship 
between ground and aerial volume. 

Discussion 

The importance of this study firstly lies on the fact that 
by implementing two-phased method using aerial pho- 
tographs and terrestrial sample plots the whole forested 



Tab. 3. Mean volume predicted by two-phased method, number of ground and aerial plots for each age class unit in the territory 
of Stolowe Mts NP 



Total no 
of aerial plots 


Reduced no 
of aerial plots 


Plots measured on the ground 
and air-photo 


Age 

classes 


Mean volume 
[m 3 ha' 1 ] 


Cor. coefficient R 


60 


51 


9 


I 


153 


0.926 


43 


37 


6 


II 


90 


0.967 


175 


154 


21 


III 


125 


0.716 


212 


184 


28 


IV 


211 


0.785 


213 


192 


21 


V 


207 


0.736 


306 


269 


37 


VI 


335 


0.775 
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R - linear correlation coefficient between the standing volume calculated on the peimanent ground control plots and on aerial 
sample plots. 
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area of Stolowe Mts NP was inventoried. The accuracy 
was similar to the result of inventory for the manage- 
ment planning (carried on ground plots only - approx. 
800 plots). Secondly, it proved that aerial photography 
could be applied practically in inventorying the moun- 
tainous forests and yield good results. Thirdly, the in- 
corporation of DTM with the aerial photography has an 
enormous effect on the overall results especially regard- 
ing to height measurements, more so when we dealing 
with mountainous forests. 

The results of this study are comparable with re- 
sults achieved in other similar previous European 
studies which used different regression equations 
in predicting the standing volume. Ranked chrono- 
logically can be mentioned results from Stellingwerf 
(1962) who had conducted the measurements over 
50-70 years old pine stand (R = 0.921); Bogyay (1970) 
over 10-60 years old pine stand (R = 0.96); Akca et 
al. (1997) achieved an (R = 0.847) over mixed stands 
of spruce and beech. The most current study to which 
the results of these study are relevant in comparison 
is Miscicki (2009), R = 0.868 whose study was con- 
ducted likewise this study over a mixed forested area. 
Needless to say the results from the other studies are 
slightly higher than the ones from the current study. 
The reason for this difference is the fact that the meas- 
urements in these studies were conducted mainly over 
stands with one dominant tree species except Miscicki 
(2009), whereas the current study comprised a forest- 
ed area with various tree species and different ages. 
Therefore the results are considered as reliable and 
comparable with the highest results found in the lit- 
erature. 

Conclusions 

- The overall results indicate that CIR aerial photog- 
raphy may be used efficiently in inventorying the 
mountainous forests. 

- Statistical analysis indicates that there were no sta- 
tistically significant differences in mean values, cal- 
culated for all age classes when using both methods 
(two-phased method and classic ground method). 

- Tree height was one of the parameters which was 
measured accurately and was identified as one of the 
key parameters in calculating the standing volume. 



- Tree height and green volume were in particular the 
most important variables in generating the final re- 
gression model, which was used in estimating the 
standing volume. 
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